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ABSTRACT 

The disruption of binary stars by the tidal field of the black hole in the 
Galactic Center can produce the hypervelocity stars observed in the halo. We 
use numerical models to simulate the full spectrum of observable velocities of 
stars ejected into the halo by this binary disruption process. Our model includes 
a range of parameters for binaries with 3-4 Mq primaries, consideration of radial 
orbits of the ejected stars through an approximate mass distribution for the 
Galaxy, and the impact of stellar lifetimes. We calculate the spectrum of ejection 
velocities and reproduce previous results for the mean ejection velocity at the 
Galactic center. The model predicts that the full population of ejected stars 
includes both the hypervelocity stars with velocities large enough to escape from 
the Galaxy and a comparable number of ejected, but bound, stars of the same 
stellar type. The predicted median speeds of the population of ejected stars 
as a function of distance in the halo are consistent with current observations. 
Combining the model with the data also shows that interesting constraints on 
the properties of binaries in the Galactic Center and on the mass distribution in 
the Galaxy can be obtained even with modest samples of ejected stars. 

Subject headings: Galaxy: center - Galaxy: halo - (stars:) binaries: general - 
stellar dynamics 

1. Introduction 

Hypervelocity stars (HVSs) are objects moving with speeds sufficient to escape the 
gravitational influence of the Galaxy. Hills (1988) predicted their existence, and Brown et 
al. (2005) discovered the first HVS in a survey of blue stars within the Galactic halo. We 
now know of seven HVSs (Edelmann et al. 2005; Hirsch et al. 2005; Brown et al. 2006a,b). 
One HVS is a subdwarf O star at a distance of ~ 20 kpc from the Galactic Center (Hirsch 
et al. 2005). The other HVSs are probably B-type main sequence stars with galactocentric 
distances of 50-100 kpc. The five late B stars in this group have masses m ~ 3-4 Mq (Brown 
et al. 2006a,b). The early B-type HVS is more massive, m ~ 8 (Edelmann et al. 2005). 

Only two proposed models plausibly explain the origin of the observed population of 
Galactic HVSs. In the single black hole model, a binary system strays too close to Sgr A*, 
the 3.5 million solar-mass black hole in the Galactic Center (Hills 1988, 1992; Yu & Tremaine 
2003; Gould & Quillen 2003). This encounter tears the binary apart, leading to the capture of 
one star and the high-speed ejection of its partner (Hills 1988). For typical ejection velocities 
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of ~ 2000 km s ^, the HVS can reach galactocentric distances of ~ 50 kpc in ~ 25 Myr 
(Hills 1988, 1992; Yu & Tremaine 2003; Gould & Quillen 2003; Brown et al. 2006b; Pcrets 
et al. 2006), roughly the lifetime of a 9 star (Schaller et al. 1992; Schaerer et al. 1993). 

Other HVS models involve three-body interactions between a single star and a binary 
black hole. Random encounters between stars near the Galactic Center and a binary black 
hole can lead to an HVS with a typical ejection velocity of ~ 2000 km (Quinlan 1996; Yu 
2002; Yu & Tremaine 2003; Gualandris et al. 2005; Scsana ct al. 2006). Alternatively, the 
in-spiral of an intermediate mass black hole with m ~ 10'^ — 10"^ Mq can eject stars from the 
Galactic Center (e.g., Hansen & Milosavljcvic 2003). Some of these encounters will produce 
HVSs with ejection velocities of ~ 2000 km s^^ (Levin et al. 2005; Baumgardt et al. 2006; 
Sesana et al. 2006). Although binary black holes with circular orbits and in-spiraling black 
holes produce roughly isotropic distributions of HVSs on the sky, eccentric binary black 
holes can produce anisotropic space distributions of HVSs (HoUey-Bockelmann et al 2005; 
Baumgardt et al. 2006; Sesana et al. 2006). 

Here, we use numerical simulations to predict the distribution of observable radial ve- 
locities of HVSs ejected by a single black hole at the Galactic Center. Our results identify a 
new population of ejected stars of the same stellar type as the HVSs but more slowly moving. 

Ongoing halo surveys (Brown ct al. 2006b) show evidence of this population. The distribu- 
tion of observable radial velocities among these stars may yield information about the nature 
of the progenitor binaries and can serve as a probe of the distribution of mass throughout 
the Galaxy. Even the small number of known HVSs can constrain the distribution of mass 
in the Milky Way (see also Gnedin ct al. 2005). 

We describe the basic model and the numerical simulations in §2, compare the numerical 
results with observations in §3, and conclude with a discussion and summary in §4. 

2. The Model 
2.1. Overview 

The goal of this study is to develop methods for predicting the observable velocity 
distribution of HVSs in the Galactic halo. We derive velocity distributions appropriate for 
the single black hole at the Galactic Center; we do not consider models for ejections from 
binary black holes (e.g., Yu & Tremaine 2003; Gualandris et al. 2005; Levin et al. 2005). 
Previous analytical and numerical analyses have concentrated on the production of HVSs 
from the interactions of binary stars with the massive black hole in the Galactic Center (Hills 
1988, 1992; Yu & Tremaine 2003; Gould & Quillen 2003; Perets et al. 2006). However, the 
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observable distribution of HVSs also depends on the mass distribution in the Galactic Center 
and on the gravitational potential of the Milky Way (e.g., Gnedin et al. 2005). Here, we 
adopt a simple prescription for the Galactic mass distribution, integrate the orbits of stars 
ejected from the Galactic Center, and derive a first approximation to the full observable 
phase-space distribution of ejected stars on radial orbits in the halo. Comparisons between 
observations of halo stars with the model predictions then yield constraints on the mass 
function of binaries in the Galactic Center and the Galactic potential. 

We simulate the ejection of stars from the Galactic Center in three steps. Following 
Hills (1988), we first integrate orbits of binary stars passing by Sgr A* to quantify ejec- 
tion probabilities and velocities. Then we use a Monte Carlo code based on semianalytical 
approximations for rapid generation of simulated catalogs of ejected stars in the Galactic 
Center. Finally we integrate the orbits of these stars in the Galactic potential to derive how 
they populate the Galaxy's halo. To generate observable samples, the simulation incorpo- 
rates the stellar main sequence lifetime from published stellar evolution calculations. In this 
first study, we focus on binaries with the full range of possible initial separations and with 
primary star masses of 3-4 M0, corresponding to halo B stars targeted in ongoing radial 
velocity surveys (Brown et al. 2006a,b). 



2.2. Simulations of Binary Interactions with the Central Black Hole 

To follow the evolution of a binary as it interacts with a massive black hole, we use a 
sixth-order, symplectic integrator (Yoshida 1990). Bromley & Kenyon (2006) describe tests 
of this code in the context of a simulation of planet formation (see also Kenyon & Bromley 
2006). For the three-body simulations here, the algorithm restricts energy errors to less than 
one part in 10^. Our numerical trials start with randomly oriented binary stars, launched 
toward the black hole from a distance of several thousand AU, much larger than the binary 
separation (generally less than a few AU). As in Hills (1988), we assume that the initial 
approach speed of the binary center of mass is 250 km s"^, and that the binary star orbit 
has negligible eccentricity. Variations in the initial approach speed have little impact on the 
results. Although the binary eccentricity can affect the outcome of an ejection event, stars in 
eccentric binaries spend most of their time at separations larger than the orbital semimajor 
axis. Thus, results for eccentric binaries generally mimic results for circular binaries with 
appropriately larger semimajor axes. Table 1 lists the ranges of other input parameters used 
in these calculations. 

Our simulations reproduce published numerical results and they are consistent with 
analytical expectations for the production of HVSs (Hills 1988, 1992; Yu & Tremaine 2003; 
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Gould & Quillen 2003). The high-speed ejection of a binary member from an encounter 
with a massive black hole depends on the binary's semimajor axis abm, the closest-approach 
distance between the binary and the black hole Rnun, and the masses of the three bodies 
— M, is the black hole mass, and mi and 1712 arc the masses of the primary and secondary 
binary partners, respectively. For encounters which produce an ejected star, our numerical 
simulations yield (see also Hills 1988) 

„, = 1760 (^r'"- (^h±^y" ( ^- km (1) 

'J VO.l AUJ V 2Mq j X lO^Moy ^ ^ 

where the factor is of order unity, and serves to tune the ejection speed in accord with 
numerical simulations (as in Fig. 3 of Hills 1988). To approximate Hills' results, we derive a 
quintic polynomial fit for /r 

fR = 0.774+ (0.0204+ (-6.23 X 10"^ + (7.62 X 10"*^ + 
(-4.24 X 10"^ + 8.62 x 10'^^ D)D)D)D)D, 

where Hills' parameter D is defined in equation (3) below. For an unequal mass binary, the 
ejection speeds of the primary and secondary are 

2m2 \ "^^^ , / 2mi ^ 



respectively. 



vi = ■ and V2 = v^j ■ , (2) 

' nil + ^2 / Vmi + 777.2 



We confirm that both stars are rarely ejected during an encounter (Hills 1988). Further- 
more, for the types of binary orbits and mass ranges we consider, the primary and secondary 
have an approximately equal chance of ejection. The key factors in the ejection are the 
orbital phase and orientation of the orbital plane when the binary comes under the grip of 
the black hole's tidal field. Our simulations suggest that if a binary is disrupted, then the 
star on the orbit closest to the black hole tends to be captured; its partner is ejected. When 
there is a significant difference between primary and secondary masses, the primary is more 
hkely to be ejected. We can understand this preference in the limit of a very large mass 
ratio, where the primary follows its unbound Keplerian orbit about the black hole and the 
secondary may get captured. However, for the mass ratios and orbital configurations which 
produce the ejection speeds of interest here, our simulations suggest that we may neglect 
any preference in the ejection of the primary star. 

The ejection speeds in equations (1) and (2) represent averages, and are the theoretical 
speeds at infinite distance from the black hole in the absence of other gravitational sources. 
Figure 1 illustrates our calculated distribution of ejection speeds for equal-mass binaries after 
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encountering Sgr A*, with M, = 3.5 x 10^ Mq. The distribution has a pronounced peak at 
~ 2800 km and broad wings extending to ~ 1000 km and ~ 4000 km s~^. The sohd 
hne in Figure 1 shows that a Gaussian with Vavg ~ 2600 km s~^ and 0"^, ~ 0.2 Vavg ~ 500 
km s~^ provides a reasonable characterization of the numerical results. Repeated trials for 
a wide range of binary configurations indicate that the Gaussian model generally gives an 
adequate description of the ejection speed distribution. We use this simple Gaussian model 
as input to the Monte Carlo simulations of the predicted distribution of speeds as a function 
of distance from the Galactic Center. Our simulations confirm that the final distribution of 
ejected stars in the halo is insensitive to our particular choice of within reasonable limits. 

The shape of the velocity distribution in Figure 1 reflects the detailed dependence of 
an ejection event on the orbital phase and orientation of the binary orbit relative to the 
black hole. We speculate that the extended tails of the distribution correspond to infrequent 
encounters where the binary angular momentum is either strongly aligned or counteraligned 
with its ccntcr-of-mass angular momentum about the black hole. Otherwise, encounters 
yield velocities near the mean, causing a peak in the distribution near Vavg- A study of a 
large ensemble of simulations may confirm this speculation and is beyond the scope of this 
initial study. 

Given values of Obin and -Rminj we derive the probability that a particular interaction 
leads to the ejection of a star. Following Hills (1988), we calculate a dimensionless quantity. 



for < D < 175. For D > 175, -Rmin ^ ^bm, and the binary docs not get close enough to 
the black hole for an ejection. Thus Pgj is zero in this case. 

Our treatment of the ejection probability does not account for the possibility that a black 
hole-binary encounter will lead to a collision between the binary partners (Yu & Tremaine 
2003; Ginsburg & Loeb 2006). Our simulations suggest that stellar collisions may be im- 
portant for small values of Obin, but only at the 0(10%) level. Furthermore, there is enough 
uncertainty in the possible outcome of such high-speed collisions, including rapidly rotating 
ejected stars or mergers (Ginsburg & Loeb 2006), that we do not consider their effect in this 
work. Instead we note that equation (4) may modestly overestimate the ejection probability 
in some cases. 

Next, we calculate the broader distribution of ejection velocities that arise from encoun- 
ters with a range of astrophysically relevant Cbin and i?min values. Surveys of large samples 




(3) 



and derive the probability of an ejection as 



Pei ~ 1 - D/175 



(4) 
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of binary systems with solar-type primary stars (e.g. Abt 1983; Heacox 1998) suggest that 
the probabihty density function for binary semimajor axis is roughly 

p(abin)rfabin ~ <iabm/abm (5) 

for abin ~ 10~^ — 10^ AU. Duquennoy & Mayor (1991) favor a Gaussian distribution in log 
P, where P is the orbital period. Here we adopt equation (5), which yields an equal number 
of objects in equally spaced logarithmic bins of semimajor axis and is consistent with the 
distribution of Obin derived from observations of smaller samples of binaries with 0-type 
and B-type primary stars (Garmany ct al. 1980; Kobulnicky et al. 2006). The probability 
density function of closest approach distances to Sgr A* varies hnearly with it^min as a result 
of gravitational focusing (Hills 1988). With these probabihties, we use a Monte Carlo code to 
draw samples of Obin and it!min- We reject some fraction of samples according to equation (4), 
and feed the rest into equations (1) and (2) to derive average ejection speeds. A normally 
distributed random number generator then provides ejection speeds with a 20% dispersion 
about the mean (eq. [1]). Thus, we produce catalogs of ejection events for specific binary 
masses, and construct distributions of ejection speeds. 

Unlike previous studies of ejected stars from the Galactic Center, which were focussed 
on HVSs, we admit wide ranges of Obin and -Rmin to consider ejection events that would not 
necessarily yield stars capable of escaping the Galaxy. For the binary semimajor axis, we 
choose 0.05 AU < abin ^ 4 AU; the lower limit is suggested by the physical radius of the 
4 M0 primary (about half this length); the upper limit arises because we are interested in 
stars with sufficient ejection speeds to populate the Galactic halo beyond 10 kpc (cf. eq. [1]). 
Once we set the maximum Obin value, equation (4) gives the maximum value of it!min that 
could result in a non-zero ejection probability. 

The semianalytical approximations in equations (l)-(5) encode the underlying physics 
of ejections from the Galactic Center. Hills' parameter D and the expression for the ejection 
probability in equation (4) are statements that the abihty of Sgr A* to disrupt a binary 
depends most strongly on the ratio between the orbital separation and the closest approach 
distance; a tight binary is hard to disrupt when it does not get close to the black hole. 
However, once a binary is disrupted, the speed of the ejected star comes mainly from the 
kinetic energy available to the star from its binary orbit, which scales as l/obin (eq. [1]). 
Finally, a broad distribution of ejection speeds arises from the interplay between binary 
separation and ejection probability. For example, a small binary separation — which can 
yield a high-speed ejection — is more likely than a large abin, yet it has a low probability 
of producing an ejection for typical values of it!min (eq. [5]). Conversely, a large binary 
separation — which leads to a low- velocity ejection — is not common; however, should a system 
with large Cbin encounter Sgr A*, then an ejection event is highly hkely. Within realistic 
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ranges of aun and Rmin, ejection speeds can vary widely. 



2.3. The Galactic Potential and the Observed Velocity Spectrum 

To construct simulated catalogs of ejected stars that are relevant to HVSs searches in 
the Galactic halo, we next consider the effect of the Galactic potential. A simple parame- 
terization of the distribution of mass in the Galaxy is 

where r is the distance to the Galactic Center, po is the central density, Oc is a "core radius," 
and the index a is around two (Binney & Tremaine 1987). Unless otherwise specified, we 
set po — 1.27 X 10^ M0 pc~^, Oc = 8 pc, and a — 2; these choices give a mass within 10 pc 
of the Galactic Center of ~ 3 x 10^ Mq, as inferred from stellar kinematics in the region of 
Sgr A* (Eckart & Genzel 1997; Schodel et al. 2003; Chez et al. 2005), and a circular rotation 
speed of 220 km s"^ in the Solar neighborhood (Hogg et al. 2005). 

To produce simulated surveys, we view the Galactic Center as a fountain of ejected stars, 
spewing out at a steady rate. We focus on primary stars of 3 M© to 4 M©, because these 
objects are the targets of ongoing B star searches. Although B stars could be secondaries, 
we assume for this work that such binaries arc relatively rare. Whether an ejected star 
becomes part of a simulated survey depends on (i) the spatial extent of the survey, (ii) the 
main sequence lifetime of the star, and (iii) the time when the star was ejected relative to 
the present. 

For the survey volume, we assign a 10 kpc inner radius and an outer radius of 90 kpc 
for 3 Mq stars and 120 kpc for 4 Mq stars, consistent with the selection in the survey of 
Brown et al. (2006a,b). To derive the time since a star was ejected, we randomly generate 
an age, T, from a uniform distribution between zero and the main-sequence lifetime tms, 
which is 350 Myr for 3 M© stars and 160 Myr for 4 M© stars (Schaller et al. 1992; Schaerer 
et al. 1993; Brown et al. 2006a). In doing so, we assume that the stars' pre-main sequence 
phase is comparatively short, that the travel time between the parent binary and the black 
hole is also short, and that the scattering of the host binary toward Sgr A* is a random 
process which is equally likely at any time (see Yu & Tremaine 2003 and Perets et al. 2006 
for descriptions of scattering mechanisms) . 

We use a second orbit integration code to calculate the radial trajectory of an ejected 
star as it travels through the potential generated by the mass density in equation (6). This 
code takes its input from the Monte Carlo ejection-speed generator, modified to account for 



-9- 



the finite distance from the black hole. We start each ejected star's orbit near Sgr A*, at 
a distance where the Galactic mass model generates 5% of the black hole mass, typically 
within a parsec. The code integrates the orbit outward from this starting point, with a simple 
leapfrog (second-order symplectic) scheme to track radius and radial velocity as a function 
of time. It concurrently integrates the density (eq. [6]) using Simpson's rule, to derive values 
of the mass enclosed at each radius, a quantity needed for force evaluation. The integrator 
breaks off when the time reaches the star's randomly chosen age, T, or if the star falls back 
out of the survey volume toward the black hole. In this way, we build up simulated catalogs 
of ejected stars within the Galactic halo. 

Figures 2-4 summarize the main results of these simulations. Figure 2 shows the observ- 
able velocity distribution as a function of secondary mass. The individual panels correspond 
to a 4 M0 primary, and secondary masses of 0.5, 1, 2 and 4 M©, respectively. In all cases, the 
velocity distribution of stars in the range 10 kpc < r < 120 kpc is peaked near 300 km s~^, 
with virtually no dependence on the mass of the secondary. Instead, the location of this 
peak is largely determined by the radial extent of the survey volume and the main-sequence 
lifetime. In the fountain picture, many of the slowly moving ejected stars can reach the inner 
halo (r < 40 kpc) within a time t < t^s, and some will even fall back toward the Galactic 
Center. However, the ejection speeds must be faster just to reach larger radii within the stel- 
lar main sequence lifetime. Thus, more distant regions of the halo tend to hold faster-moving 
ejected stars. Figures 2-4 show this effect by comparing the velocity distribution peak for 
stars inside r = 40 kpc (blue histograms) with the distribution for stars located outside of 
that radius (magenta histograms). The velocity distributions for nearby stars peak at ~ 
100 km s""^; the peak moves to 400-500 km for more distant objects. 

The overall velocity distribution of primary stars is insensitive to the mass of the sec- 
ondary. In general, the stellar hfetime combined with the Galactic potential 'filter' the 
distribution to admit slower-speed ejected stars which are prevalent independent of the sec- 
ondary mass. These slower-speed ejecta live long enough to populate the inner Galactic 
halo in larger numbers compared to the high velocity tail. At larger distances, only the high 
velocity tail is observable as a result of the short stellar lifetimes. The dependence of ejection 
velocity on secondary mass (equations 1-2) is important only for the high velocity tail of the 
observed velocity distribution. For example, the observation of a blue star with a speed in 
excess of 2000 km s~^ requires a secondary of at least 1 M©, and statistically favors an even 
larger mass because the likelihood of a hypervelocity ejection increases with the secondary 
mass. 

Figure 3 illustrates the effect of the Galactic mass density profile on speeds of ejected 4 
M0 stars in the halo. Although the peak of the velocity distribution is relatively independent 
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of the mass model, the median velocity is sensitive to model parameters. Changing the index 
a in equation (6) from 1.8 to 2.5 produces a significant variation in the median velocity of 
objects in the full survey volume (black histograms in each panel), from 270 km s~^ for a 
= 1.8 to more than 600 km s~^ for a = 2.5. The blue and magenta histograms show the 
sensitivity of the median to distance from the Galactic Center. For halo stars beyond 40 
kpc, the median velocity varies from 485 km s~^ for a = 1.8 to 750 km s^^ for a = 2.5. 

Figure 4 shows similar results for binaries with a primary mass of 3 Mq and secondary 
mass of 1.5 Mq. Compared to the 4 Mq cases, the histograms shift to slower ejection speeds 
due to the lower luminosity and longer (350 Myr) lifetime of the lower mass stars. At 
the magnitude limit of optical surveys, stars with smaller optical luminosities have smaller 
distances; thus these surveys are shallower for 3 Mq stars than for 4 Mq stars. 

Table 2 summarizes the results for the median radial velocities in Figures 3 and 4. For 
each adopted a, the Table lists the core radius (ac) and central mass density (po) needed for 
a mass of ~ 3 x 10^ Mq within 10 pc of the Galactic Center and the median velocities for 
two samples of ejected stars for 3 Mq and 4 Mq primary stars. For 3 Mq primaries, the full 
survey volume extends from 10 kpc to 90 kpc; the halo volume extends from 30 kpc to 90 
kpc. For 4 Mq primaries, the full survey volume extends from 10 kpc to 120 kpc; the halo 
volume extends from 40 kpc to 120 kpc. 

In addition to producing smaller median observable speeds for lower mass stars, the 
main sequence lifetimes set the maximum stellar mass for HVSs. For primary stars ejected 
with velocity Vgj from the Galactic Center, Table 3 lists the maximum distance dmax the 
star can reach during its main sequence lifetime. Because the travel time from the Galactic 
Center to the Galactic halo, t oc m~^/^, has a much weaker mass dependence than the main 
sequence lifetime, tms m^^, massive stars ejected from the Galactic Center cannot reach 
large galactocentric distances. Furthermore, high ejection speeds require a small binary 
semimajor axis (eq. [1]); larger, more massive stars have a larger minimum abin than smaller, 
less massive stars, based on an assessment of Roche lobe overflow (Eggleton 1983). Thus, 
less massive stars have comparable maximum ejection velocities - and hence larger d^ax, 
given enhanced stellar hfetimes - than more massive stars. Based on the relative ejection 
velocities, stellar lifetimes, and the time needed for any star to reach the central black hole, 
we suggest that the maximum mass for an HVS in the outer halo is ~ 12 Mq (Table 3). 
With large surveys to greater depth, this mass limit may provide an observational test of 
these models. 

Figures 2-4 show the breadth of speeds that ejected stars can have once they have 
traveled into the Galactic halo, from negative values — corresponding to infall — to hy- 
pervelocities. Clearly, HVSs have slower-moving counterparts of the same stellar type that 
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originated near the Galactic Center and carry information about the physical process which 
led to their ejection. To illustrate this feature of the simulations in more detail, Figure 5 
shows the fraction of escaping stars (with velocities exceeding the local escape velocity^) as 
a function of distance from the Galactic Center. For 3 Mq primaries, the fraction increases 
roughly linearly with radius and is almost 50% at ~ 90 kpc, the limit of current surveys for 
these stars. The fraction of escaping stars is much larger for 4 Mq primaries and is 100% 
at 80-120 kpc. In an ensemble of 3-4 Mq ejected stars, the fraction of escaping stars is ~ 
50%, suggesting that observations should reveal roughly comparable numbers of unbound 
HVSs and ejected, but bound, halo stars with lower observed radial velocities. To test this 
and other aspects of the calculation, we now compare our results with observed velocity 
distributions of stars in the Galactic halo. 



3. Comparing the model with observations 

The model for HVSs has several observable and testable consequences. To produce an 
observable ensemble of ejected stars, the Galactic Center must have a significant population 
of young binary stars close enough to the central black hole. 

Current observations provide broad support for many young stars at the Galactic Center 
(e.g., Martins et al. 2005). Near Sgr A*, high spatial resolution observations reveal ~ 100 
or more young OB stars and many evolved Wolf-Rayet stars (Chez et al. 2003; Schodel et 
al. 2003; Eiscnhauer et al. 2005; Tanner et al. 2005; Paimiard et al. 2006). Within 25- 
50 pc of Sgr A*, there are several distinct star-forming regions, including the Arches cluster, 
containing thousands of stars with masses > 3-4 Mq (Figer et al. 1999; Najarro et al. 2004; 
Stoltc et al. 2005). If the binary fraction of the OB population near the Galactic Center 
is comparable to the local fraction of ~ 70% (Kobulnicky et al. 2006), then the Galactic 
Center contains enough binary stars to interact with the central black hole and to produce 
an observable population of HVSs. 

The ejection rate from the central black hole depends on the structure of the inner 
Galaxy near Sgr A*. For a random phase-space distribution of binaries in the Galactic 
Center, Hills (1988) first argued that the ejection rate could be as high as ~ 10~^ yr~^- Yu 
& Tremaine (2003) refined the calculation to account for the destruction of binaries whose 
orbits take them near Sgr A*. This region of phase-space — the loss cone — is replenished by 
random scattering between binaries and other stars. If a steady state is reached, the ejection 



^We calculate the escape velocity from the mass distribution in equation (6) for a Galaxy with an outer 
radius of 120 kpc. The results are nearly identical for outer radii of 250 kpc. 
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rate is ~ 10~^ yr~^ for objects with binary separations within 0.3 AU. Recently, Perets et 
al. (2006) conclude that phase-space mixing by massive star clusters, molecular clouds, and 
perhaps intermediate mass black holes produces a significantly larger ejection rate. They 
estimate a Galactic population of 10-200 HVSs with stellar masses of 4 Mq, consistent with 
observations (Brown et al. 2005, 2006a). 

Once these ejected stars arc in the halo, the velocity distributions derived from ongoing 
radial velocity surveys provide a test of the ejection model. To make a first comparison, we 
use data from the halo survey of Brown et al. (2006a,b) and data for 2 HVSs from Edclmann 
ct al. (2005) and Hirsch ct al. (2005). The Brown ct al. survey selects candidate B-type 
stars from the Sloan Digital Sky Survey (SDSS; Adchnan-McCarthy et al. 2006) and derives 
distances, radial velocities, and spectral types from optical spectra acquired with the 6.5-m 
MMT. This survey is complete and can detect all B-type stars at any velocity out to ~ 100- 
150 kpc. In addition to HVSs with radial velocities exceeding ~ 400 km s~^. Brown et al. 
(2006b) find an overabundance of B-type halo stars with radial velocities of ~ 250-400 km 
s~^. Because the observations show no evidence for an infalling population with comparable 
radial velocities. Brown et al. conclude that these outliers in the velocity distribution might 
be low velocity HVSs or the high velocity tail of runaway B-type stars (Portegies Zwart 
2000). 

To compare our predictions with the observations. Figure 6 shows the median radial 
velocity as a function of radial position in an a = 2.0 model for 3 Mq and 4 Mq primaries 
with secondary star masses equal to half the primary mass. The radial variation in the median 
speed refiects the interdependence of the ejection velocity distribution, stellar lifetime, and 
survey volume. The trend toward higher median speed at larger results from the finite stellar 
lifetime of ejected stars; slow-moving stars cannot reach large distances within their main- 
sequence hfetimes (see also Baumgardt et al. 2006). With a sizable population of ejected 
star candidates, the velocity-radius relationship in this Figure serves as a fundamental test 
of the stellar ejection hypothesis. 

The comparison of the model prediction with observations of known high velocity halo 
stars (open symbols in Figure 6) is encouraging. In general, the model predicts roughly 
comparable numbers of HVSs and stars with radial velocities of 250-400 km s~^, as observed 
in current surveys. Aside from the single subdwarf O star outlier (Hirsch et al. 2005) at 
(r, Vmed) ~ (20 kpc, 700 kms), the observed velocity distribution of HVSs also closely follows 
the predicted relation for 4 Mq primary stars. The lone 8 Mq HVSs (Edelmann et al. 2005) 
at (r, Vmed) ~ (55 kpc, 550 kms) lies within the group of lower mass B-type stars^. At 



^We include the sdO star and the 8 M© B star for completeness. Our predicted distribution does not 
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lower velocities, the observed distribution lies between the predicted relations for 3 Mq and 
4 Mq stars. With measured spectral types consistent with stellar masses of 3-4 Mq for most 
HVSs and most of the low velocity sample, we conclude that the observations are in excellent 
agreement with model predictions. 

Although Galactic Center ejections are necessary to explain the observed radial velocities 
above 400 km s~^, objects with radial velocities of 250-400 km s~^ could be runaway B 
stars from supernova-disrupted binaries or other stellar dynamical interactions (e.g., Blaauw 
1961; Poveda et al. 1967). However, for radial velocities exceeding ~ 200 km s^^, the 
velocity distribution of runaway stars should probably fall off more steeply than is observed 
(Leonard 1991; Portegies Zwart 2000). Current model predictions and statistics of halo stars 
with radial velocities of 250-400 km do not provide a conclusive test of the runaway B 
star scenario. A simulation of the observable population of runaway B stars — similar in 
spirit to our simulation of the observable population of HVSs — would provide a good test 
of this mechanism. 

Despite the small number of plausibly ejected stars in the sample of Brown et al. (2006b), 
these objects provide an interesting test of models of the Galactic mass distribution. We 
consider the model with the steepest large-radius fall-off in density, corresponding to a = 2.5 
in Table 1, and the HVSs listed by Brown et al. (2006b, Table 1 therein), and shown here 
in Figure 6. None of the nine late-type B stars with distances in excess of 40 kpc from 
the Galactic Center and speeds above 250 km s~^ have a speed in excess of the model's 
predicted median velocity of 752 km s~^. Thus, we can rule out the model at a confidence 
level of 99.9%. If the objects with speeds less than 400 km s~^ are not ejected stars, then the 
confidence level drops to about 97%. Including the more massive early- type HVS (Edelmann 
et al. 2005), which should have a higher ejection speed than the late B-type stars, raises the 
confidence level above 98%. A similar argument using all ten B stars suggests that the 
a — 2.2 model is excluded with better than 98% confidence. The models with shallower 
large-radius fall-off can not be ruled out in this way. 

These constraints on Galactic mass models depends heavily on the assumption that 
the known late-type HVSs have masses close to 4 Mq. If all objects were 3 Mq stars, the 
predicted median velocity for ejected stars outside of 30 kpc is 437 km s~^ in the a = 2.5 
model (Table 2). This median speed is consistent with the data. However, the model predicts 
that 25% of stars should have a speed greater than about 750 km s~^. All eleven late-type 
B stars shown in Figure 6 have speeds in excess of 250 km s~^ and distances greater than 
30 kpc, yet none has a speed exceeding 750 km s~^. In this case, the a = 2.5 model may be 



include massive B main sequence stars or sdO stars. 
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ruled out at the 95% confidence level; this level is even higher if we consider the early-type 
HSV as well. Note that these results depend mainly on the mass of the ejected primary, not 
on the mass of the secondary star (see Fig. 2). 

In making these first probes of the Galactic mass distribution with HVSs, we note that 
larger samples of halo stars, globular clusters, and dwarf galaxies currently place better 
constraints on the mass distribution in the Galactic halo (e.g., Battaglia et al. 2005; Dehnen 
et al. 2006, and references therein). In particular, these data rule out the a > 2.2 models 
with higher significance than the known HVSs and lower velocity ejected stars. However, 
our results show that even small samples of HVSs are a potentially important probe of the 
Galactic mass distribution to large radii. 

4. Discussion and summary 

We have developed a first calculation of the predicted distribution of radial velocities 
for HVSs at galactocentric distances of 10-120 kpc in the Galactic halo. Our simulations 
quantify the ejection probabilities and velocities of binaries interacting with Sgr A*, generate 
the distribution of ejected stars from the vicinity of the Galactic Center, and produce the full 
phase-space radial velocity distribution of ejected stars in the halo from an integration of the 
orbits of ejected stars through the Galactic potential. In addition to providing a framework 
for interpreting the observed velocity distribution of HVSs and related stars, the results 
of these simulations (Figures 1-4) demonstrate that this approach can provide interesting 
constraints on the Galactic potential and on the mass function of the population of binary 
stars in the Galactic Center. 

Besides confirming Hills' (1988, 1992) results for the average velocity of ejected stars 
from the black hole, we derive a broad distribution of ejection velocities for each value of Obin 
and i?inin- A Gaussian with ^ 0.2vavg provides a reasonable approximation to the derived 
velocity distribution. We plan to consider the origin of this distribution in future studies. 

The initial comparisons between the model and observations are encouraging. The 
measured radial velocities of HVSs and stars with somewhat lower speeds agree well with 
the predictions. The current population of these stars provides a constraint on the Galactic 
mass distribution. As ongoing radial velocity surveys yield more HVSs with accurate stellar 
masses, the model will allow more robust tests for the origin of HVSs and lower velocity 
stars, will provide limits on the initial mass function of primary stars in the Galactic Center, 
and will yield better constraints on mass models for the Galaxy. 

As surveys reveal more HVS, the observed velocity distribution of ejected stars in the 
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Galactic halo may help to distinguish among alternative scenarios for producing HVSs. Com- 
pared to the single black hole models, binary black hole models produce relatively fewer 
ejections at the highest velocities (Yu & Tremaine 2003; Baumgardt et al. 2006). Although 
binary black hole models may yield larger ejection rates overall (e.g., Yu & Tremaine 2003), 
the two models produce similar velocity distributions at lower velocities. Thus, better con- 
straints on the HVS distribution at the highest velocities may provide a test of the two 
pictures for HVS formation. 

Refinements to the mass models can include more detailed Galactic structure, in par- 
ticular the Galactic disk. Indeed, Gnedin et al. (2005) have already proposed that HVSs can 
serve to probe halo triaxiality. Including the disk potential should lead to nonradial ejection 
velocities, and may cause some clustering of ejected stars in the halo. We may therefore use 
the breakdown of spherical symmetry in the distribution of ejected stars about the Galactic 
Center to reveal the Galaxy's structure. 

If ejected stars are produced in short-lived star forming regions near Sgr A*, or if they 
originate in individual star clusters which settle dynamically into the Galactic Center, ejected 
stars — including HVSs — may also cluster in time. Given a Galactic mass model, we 
can identify any temporal correlations by using the radial position and velocity of stars to 
determine the time since they were ejected. The extraction of temporal and spatial clustering 
information may be difficult with rare HVSs. However, a clustering analysis might be feasible 
with the more numerous, slower population of ejected stars suggested by our model. 

We thank an anonymous referee for raising interesting issues that improved the paper. 
We acknowledge support from the NASA Astrophysics Theory Program through grant NAG5- 
13278. 
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Table 1. Model parameters 



parauK'tor 


H\aul)()] 


\i\hiv or raii,t!,v 


Sgr A* mass 


M. 


3.5x10*^ Mq 


primary mass 


mi 


3-4 M0 


secondary mass 


m2 


0.5-4 M0 


binary semimajor axis 


«bin 


0.05-4 AU 


closest approach distance 


-^min 


~l-700 AU 


Galactic density index^ 


a 


1.8-2.5 


Galactic "core" radius 




1-33 AU 


Galactic central density^ 


Po 


6270-162000 Mg/pc^ 



^The parameters a, Uc and po are defined in eq. (6) 
''These values do not include the mass of Sgr A* 



Table 2. Median observable speeds of ejected stars as a function of mass model 



3 Mra'' 4 Mp,^ 



a'' 


Oc (pc) 


Po (Mo/pc^) 


Vmed (full) 


Vmed (fialo) 


Vmed (full) 


Vmed (fialo) 


1.8 


1 


1.62 X 10^ 


144 


188 


270 


485 


2.0 


8 


1.27 X 10* 


218 


272 


374 


565 


2.2 


16 


8920 


306 


370 


484 


648 


2.5 


33 


6270 


407 


457 


608 


752 



^These columns list results for a primary mass of 3 Mq and a secondary mass of 1.5 Mq. 
The radial extent of the full survey - Vmed (full) - is 10-90 kpc; the radial extent of the 
halo survey - Vmed (halo) - is 30-90 kpc. The listed values correspond to the median 
radial velocity of the ensemble of ejected stars. 

''These columns list results for a primary mass of 4 Mq and a secondary mass of 2 Mq . 
The radial extent of the full survey - Vmed (full) - is 10-120 kpc; the radial extent of the 
halo survey - Vmed (halo) - is 40-120 kpc. The listed values correspond to the median 
radial velocity of the ensemble of ejected stars. 

'^a, ttc and po are defined in eq. (6). 
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Table 3. Maximum distance from the Galactic Center as a function of ejected mass 



mi [Mq) 


J. /A JT ^\ 

tms (Myr) 


/ TV /T \ 

m2 [Mq) 


abin (AUj 


Vej (km s ) 


dmax (kpc) 


7 


43 


2 


0.04 


4500 


193 


7 


43 


4 


0.05 


4400 


184 


7 


43 


7 


0.05 


4800 


200 


9 


26 


2 


0.05 


4400 


114 


9 


26 


4 


0.06 


4200 


110 


9 


26 


9 


0.06 


4700 


123 


12 


16 


3 


0.06 


4400 


71 


12 


16 


6 


0.06 


4700 


76 


12 


16 


12 


0.06 


5200 


83 


15 


12 


4 


0.06 


4800 


58 


15 


12 


8 


0.07 


4700 


57 


15 


12 


15 


0.07 


5200 


62 


20 


8 


5 


0.07 


4800 


39 


20 


8 


10 


0.08 


4800 


39 


20 


8 


20 


0.09 


5035 


41 


30 


5 


8 


0.10 


4700 


23 


30 


5 


15 


0.11 


4700 


24 


30 


5 


30 


0.12 


5000 


25 



'^These columns list the average ejection velocity Vej for the specified binary 
semimajor axis abin, and the corresponding maximum galactocentric distance 
dmax that an ejected primary can achieve, as a function of the primary mass mi 
and the secondary mass m2. Specifically, if the primary is ejected with velocity 
Vej, dmax is thc distaiicc from the Galactic Center the star reaches during its 
main sequence lifetime fms as it travels radially outward through thc Galactic 
halo. Here, Vej is taken from eq. 1 with = 1, that is, assuming an optimal 
impact parameter -Rmin- The binary semimajor axis abin is chosen to reflect a 
minimum separation as a result of Rochc-Iobc overflow. Note that it would be 
rare for a primary to actually achieve the listed dmax^ 

since ejections are not 

common for very tight binaries. 
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Fig. 1. — The distribution of ejection speeds from a binary with 4 M0 primary stars, just 
after encountering a black hole with a mass of 3.5 x 10® Mq. The original binary has a 
semimajor axis of 0.1 AU, and its center of mass is targeted to reach a minimum distance 
of 5 AU from the black hole. 
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Fig. 2. — The distribution of observable speeds of ejected 4 primaries in the Galactic 
halo, as a function of the mass of the secondary. The secondary mass is indicated in each 
panel. The solid black histograms show the velocity distribution of all objects in a survey 
volume which extends from 10 kpc to 80 kpc away from the Galactic Center. The gray 
vertical lines locate the median velocity of each distribution. The blue histogram represents 
all objects inside of a 40 kpc radius; the magenta histogram corresponds to objects outside 
this radius. 
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Fig. 3. — Same as the previous figure, but for tlie distribution of observable speeds of an 
ejected 4 Mq primary with a 2 M0 secondary, as a function of Galactic mass model. The 
model parameters, displayed in each panel, correspond to eq. (6). 
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Fig. 4. — Similar to the previous figure, but for tlie distribution of observable speeds of 
an ejected 3 Mq primary with a 1.5 Mq secondary, as a function of Galactic mass model. 
The survey volume in this case spans a radial extent of 10 kpc < r < 90 kpc; the blue and 
magenta histograms correspond to populations inside and outside of 30 kpc, respectively. 
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Fig. 5. — The fraction of observable ejected stars with radial velocities exceeding the escape 
velocity as a function of radial distance from the Galactic Center. As indicated in the legend, 
squares show results for 4 M© primary stars; triangles show results for 3 Mq primaries. 
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Fig. 6. — The median speeds of ejected 3 (filled triangles) and 4 M0 (filled squares) 
primaries in the Galactic halo, as a function of distance from the Galactic Center. In 
both cases the secondary mass is half that of the primary. The open boxes are the five 
HVSs discovered in a targeted survey (Brown et al. 2006b); the open triangles are lower 
velocity, plausibly ejected, stars from this survey. The open circles are 2 HVSs discovered 
independently (Edelmann et al. 2005; Hirsch et al. 2005) and which are outside the mass 
range we consider. 



